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ABSTRACT 

We develop a useful formula for power spectrum analysis for high and inter- 
mediate redshift galaxy samples, as an extension of the work by Feldman, Kaiser 
& Peacock (1994). An optimal weight factor, which minimizes the errors of the 
power spectrum estimator, is obtained so that the light-cone effect and redshift- 
space distortions are incorporated. Using this formula, we assess the feasibility 
of the power spectrum analysis with the luminous red galaxy (LRG) sample in 
the Sloan Digital Sky Survey as a probe of the equation of state of the dark en- 
ergy. Fisher matrix analysis shows that the LRG sample can be sensitive to the 
equation of state around redshift z=0.13. It is also demonstrated that the LRG 
sample can constrain the equation of state with (1-sigma) error of 10% level, if 
other fundamental cosmological parameters are well determined independently. 
For the useful constraint, we point out the importance of modeling the bias 
taking the luminosity dependence into account. We also discuss the optimized 
strategy to constrain the equation of state using power spectrum analysis. For a 
sample with fixed total number of objects, it is most advantageous to have the 
sample with the mean number density 10~ 4 /i 3 Mpc~ 3 in the range of the redshift 
0.4 < z £ 1. 

Subject headings: cosmology : theory — galaxies: clusters: general — large scale 
structure of universe 
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1. INTRODUCTION 

The clustering of the cosmological objects like galaxies, clusters of galaxies and QSOs, 
is the fundamental to probe the Universe because it directly reflects properties of dark 
components and the primordial density fluctuations. The power spectrum is a simple but 
very useful tool to characterize their spatial distribution. Actually, useful constraints on the 
cosmological parameters are obtained from the power spectrum analyses of the Two Degree 
Field (2df) Galaxy Redshift Survey (Percival et al. 2001) and the 2df QSO Redshift Survey 
(Hoyle et al. 2002). In such power spectrum analyses, many authors base their methods on 
the seminal paper by Feldman, Kaiser & Peacock (1994, Hereafter FKP). 

For redshift surveys such as the 2df survey and the Sloan Digital Sky Survey (SDSS), 
however, several observational effects on the power spectrum analysis can be very important. 
Because cosmological observations are feasible only on the light-cone hypersurface defined 
by the current observer, the effect of the redshift evolution of the luminosity function, the 
clustering amplitude, and the bias, contaminates the observational data. We call this the 
light-cone effect (Matarrese et al. 1997; Matsubara, Suto & Szapudi 1997; de Laix & Stark- 
man 1998; Yamamoto & Suto 1999). On the other hand it is well known that the distribution 
of the sources in redshift sapce is different from that in the real space due to the redshift- 
space distortions. The linear redshift distortion is the effect of the bulk motion of the sources 
within the linear theory of density perturbation (Kaiser 1987; Hamilton 1998 and references 
therein). The finger of Got effect is the redshift distortion due to the random motion of the 
sources in the nonlinear regime (Mo, Jing & Boerner 1997; Magira, Jing & Suto 2000). The 
geometric distortion is the effect caused by a choice of the distance-redshift relation to plot 
a map of the sources (Ballinger, Peacock & Heavens 1996; Matsubara & Suto 1996). 

The LRG sample of the SDSS spectroscopic survey will provide a sample of 10 5 intrin- 
sically luminous early-type galaxies to z ~ 0.5 (Eisenstein et al. 2001). In the analysis of 
clustering statistics of such a sample, the light-cone effect and the redshift-space distortions 
can be substantial. Fortunately these observational effects have been well investigated and 
we can model the power spectrum incorporating them (see section 3.1). The primary pur- 
pose of the present paper is to extend the formulas in FKP in order to incorporate these 
observational effects and to derive a generalized expression of the optimal weight factor in 
the power spectrum analysis (cf. Tegmark 1995). 

The results of the 2df galaxy survey (Peacock et al. 2001), weak lensing surveys (Re- 
fregier et al. 2002; Bacon et al. 2002; and references therein), cosmic microwave background 
anisotropy measurements (e.g., Ruhl et al. 2002) and type la supernovae measurements 
(Perlmutter et al. 1999; Riess et al. 1998), support the concordance model (Wang et al. 
2000): A spatially flat universe dominated by the dark energy component and, with respect 
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to structure formation, cold dark matter with the primordial density fluctuations predicted 
in the inflationary scenario. The very recent result by the Wilkinson Microwave Anisotropy 
Probe (WMAP) strongly supports the concordance model (Spergel et al. 2003). Now the 
mystery of the dark energy component has become one of the most important issues in cos- 
mology, which has led to recent activity in investigating the dark energy (See e.g. Peebles 
& Ratra 2003, for a recent review, and references therein). The dark energy can be charac- 
terized by its equation of state wx = Px/px, where px is the pressure and px is the energy 
density. For the cosmological constant, w x = — 1- However, if the dark energy originates 
from the vacuum energy of a variable scalar field, like the quintessence model, wx can take 
values wx > — 1, and can in general be a function of redshift. Thus constraints on the equa- 
tion of state is quite important in considering the origin of the dark energy. Then various 
strategies for probing the equation of state have been investigated (e.g. Newman & Davis 
2000, Saini, et al. 2000, Wang, et al. 2000, Chiba & Nakamura 2000, Huterer & Turner 2001 
, Yamamoto & Nishioka 2001, Kujat et al. 2002, and references therein). 

Second purpose of this paper is to assess the feasibility of measuring the equation of 
state using the power spectrum analysis. 1 Recently, Matsubara & Szalay have discussed the 
usefulness of the LRG sample in SDSS for measuring the cosmological parameters (2002). 
Their method is based on maximum likelyhood analysis in redshift space. They demonstrated 
the usefulness of various SDSS samples to constrain cosmological parameters by estimating 
the Fisher matrix element. Motivated by their work, we assess the feasibility of the method 
with the power spectrum analysis of the LRG sample in SDSS to probe the equation of state 
based on the Fisher matrix formalism for power spectrum analysis (Tegmark 1997, Tegmark 
et al. 1998a). The advantage of our method is that the formulae are simple and analytic, 
which allows us to evaluate the formulae easily and to understand meaning of results clearly. 

This paper is organized as follows: In section 2, we derive formulae for the power 
spectrum analysis taking the light-cone effect and the redshift-space distortions into account. 
In section 3, the Fisher matrix element is evaluated using an approximate method. There 
we focus on the matrix elements relevant to a measurement of dark energy, especially the 
equation of state. In section 4, we discuss optimizing the strategy, for a sample with a fixed 
total number of objects, using the power spectrum to probe the equation of state. Section 5 
is devoted to summary and conclusions. Throughout this paper we use a system of units in 
which the velocity of light c equals 1. 



After completing this work, a similar investigation by Matsubara & Szalay (2003) has been announced. 
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2. FORMALISM 

In this section we present an optimal weighting scheme of the power spectrum analysis 
for redshift surveys taking the redshift-space distortions and the light-cone effect into ac- 
count. We obtain the generalized optimal weight factor, which minimizes errors of a power 
spectrum estimator, as a simple extension of the scheme developed by FKP. On the basis of 
this result, we derive the expression of the Fisher matrix element in subsection 2.2. 

In a redshift survey, the redshift z is the indicator of the distance, therefore, we need to 
assume a distance-redshift relation to plot a map of objects. For this distance, in the present 
paper, we adopt the following expression, which is the comoving distance in spatially flat 
universes 

s(z) = j_ r dz> a) 

Ho Jo v/fi (l + ^) 3 + i-^o' 
where H = 100/ikm/s/Mpc is the Hubble parameter and we fix Q = 0.3. Our formula 
presented below is general, and it does not depend on this choice of s(z). 



2.1. Optimal Weight Factor 

The derivation of the optimal weight factor in this subsection is essentially same as that 
of FKP. Some parts of equations here can be found in FKP, which we have not omitted for 
being self-contained of the present paper. We denote the number density field of sources (real 
catalog) by n g (s), where s(= s(z)^f) is the three dimensional coordinate in the (cosmological) 
redshift space and 7 is the unit directional vector. We use n(s) to denote the expected mean 
number density. We define the fluctuation field to be 

F(s) = n g (s) - an s (s), (2) 

where n g (s) = Yli^( s ~ s «)> with Sj being the location of the %— th object, similarly n s (s) is 
the density of a synthetic catalog which has mean number density 1/a times that of the real 
catalog. For n g (s) and n s (s), following FKP, we assume 

(n g (s 1 )n g (s 2 )} = n(si)n(s 2 )(l + f (si, s 2 )) + n(si)<J(si - s 2 ), (3) 

(ra a (si)n s (s 2 )) = aT 2 n(si)n(s 2 ) + a' l n(s 1 )5(s 1 - s 2 ), (4) 

(n g (s 1 )n s (s 2 )) = a _1 n(si)n(s 2 ), (5) 

where £(si,s 2 ) denotes the correlation function. 

We define the Fourier coefficient of -F(s), with a conventional weight factor, to be 

W [/rfsn(s) 2 W (s,k) 2 ]V2' [ -° ) 
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where w(s, k) is the weight function which can be adjusted to optimize the power spectrum 
estimator below. Then the expectation value of the square |jF(k)| 2 is written, 

/ dgi / ds2w(s u k)w(s 2 , k) (F( Sl )F(s 2 ))e- k -^-^) 
{l Wl ; Jcisn(s) 2 W (s,k) 2 

(7) 

Using equations (3)-(5), we have 

(F(si)F(s 2 )> = n(si)n(s 2 )e(si, s 2 ) + (1 + a)n(s 1 )5(s 1 - s 2 ). (8) 

Under the distant observer approximation, |si — s 2 | <C |si|, |s 2 |, the correlation function can 
be expressed as 

e(Sl ' S2) = / (03 P (k>l«D^ (si - S2) , (9) 

where P(k, |s|) is the redshift-space power spectrum at the distance |s|, where we defined 
s = (si + s 2 )/2. We also adopt the approximation 

w(si, k)w(s 2 , k)n(si)n(s 2 ) ~ w(s, k) 2 n(s) 2 , (10) 

which yields 

/l-ryi^|2\ _ J dsn{s) 2 w(s,k) 2 P{k, \s\) + (1 + a) J dsn{s)w{s,k) 2 

J dsn(s) 2 w(s,k) 2 • {U) 

The last term of the right hand side of equation (11) corresponds to the shot-noise, therefore 
the estimator of the power spectrum should be defined as 

P(k) = |^(k)| 2 -P shot (k), (12) 

where P s hot is defined 

_ (1 + a) J dsn(s)w(s,k) 2 

The angular average of V{k) over a thin shell of the radius k(— |k|) in the Fourier space 
gives the estimator of the angular averaged power spectrum 



Mk) = ±r I dkV(k), (14) 
>v k 



V k 
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where 14 denotes the volume of the shell. Thus the expectation value of Vo(k) gives the 
angular averaged power spectrum incorporating the light-cone effect 



dsn(s) 2 w(s,k) 2 P (k, |s|) 
(Vo(k)) = J - , , (15) 



J dsn(s) 2 w(s,kf 



P (k, |S|) = ±r I dkP(k,| S |). (16) 



where 



V kJVk 

Next we consider the variance of V(k), defined by 

(AV(k)AV(k')) = ([P(k) - (V(k))][P(k') - (V(k'))\) 

= (V(k)V(k'))-(V(k))(V(k')), (17) 

which is rephrased, by using (12), 

(AP(k)AP(k')) = (|.F(k)| 2 |.F(k')| 2 > - (|^(k)| 2 )(|^(k')| 2 ). (18) 
Substituting (6) into equation (18), we have 
(AP(k)AP(k')) = 

mlAl dsM^MUUW ^^(s J ,kQ](F(s 1 )F(s 3 ))(F(s 2 )F(s 4 ))e- k -^^)e- k/ -^^) 

[J dsn{s) 2 w{s, k) 2 ] [/ ds'n(s') 2 w(s', k') 2 ] 

(19) 

where we assumed the following relation, which is exact when F(s) follows the statistics of 
the Gaussian random field, 

(F( Sl )F(s 2 )F(s 3 )F(s 4 )) = (F(si)F(s 2 ))(F(s 3 )F(s 4 )> 

+ (F( Sl )F(s 3 ))(F(s 2 )F(s 4 )) 
+ (F( Sl )F(s 4 ))(F(s 3 )F(s 2 )). (20) 

In a similar way to derive equation (11), using the distant observer approximation, one can 
have the relation 

dsi J rfs 3W (s 1 ,k)«;(s3,k / )(F(s 1 )F(s 3 ))e tk - Sl - lk '- S3 

~ [ dsn(s) 2 u>(s, k)w(s, k')e rs < k ^ \p (^-^-, |s|) + ^1] . (21) 
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Then, with a repeated use of the distant observer approximation, we have 

2 / dsn(s) 4 w(s, k) 4 (P(k, s) + (1 + a)/n(s)f (2n) 3 5^(k - k') 



(AP(k)AP(k')) ~ 



[ f dsn(s) 2 w(s, k) 



212 



(22) 



From (14), the variance of Vo(k) is obtained by 

(APo(A:) 2 ) = (p>„(*) " (Vo(k)}] 2 ) 

= / dk f dk'(AP(k)AP(k')>, (23) 



which reduces to 



where we defined 



/AV (kV\ - 9 (2?r)3 /^(s) 4 ^(s,k) 4 g 2 (fc,s) 



Q 2 (k,s) = ^[ dk[p(k,|s|) + ^±If. (25) 
Vfc L V / n(s) J 

The optimal weight factor w(s, k) which minimizes (AVo(k) 2 ) is found from the stationary 
solution against the variation w — » w + Sw, to be 

w(s,k)= 1 (26) 
n(s)Q(k,s) 

This generalized formula incorporates the redshift distortions. In general, the redshift space 
power spectrum P(k, s) is written, in terms of the multipole expansion, 

P(k,s)= Pi(k,s)£i(»), (27) 

1=0,2,- 

where \i is the directional cosine between the line of sight 7 and the wave number vector, 
pi = cos(7 • k/fc), Ci(x) is the Legendre polynomial of the /-th order, and Pi(k, s) are the 
expansion coefficients. Then (25) reduces to 

1 2 

Q 2 (k, g) = (P (fc, |s|) + — ) + P 2 (fc, |s|) 2 + P 4 (fc, |s|) 2 + • • • , (28) 

where we considered the limit a < 1. In practice, the contribution of higher moments / > 2 is 
not large (cf. Yamamoto 2003). Neglecting this higher moments and the redshift dependence 
of the power spectrum, (26) reduces to the similar expression by FKP (cf. Peebles 1980), 

w ( s , k ) - 7 } n n 1 1 \ ? (29) 

v ; l + n(s)P (/c, |s|)' v ; 
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then, from equation (24), we have 

(AV (k) 2 } = 2^- \ [ ds- , ^ - 1 ~\ (30) 

V K J 1 V k lj (l+n(s)P (fc,|s|)) 2 J 1 ; 

Here we mention a technical problem of the above optimal weight factor. Namely, this 
weight factor contains the power spectrum at each redshift, i.e. Po(k,s[z]). It might be 
difficult to evaluate it with a sufficient accuracy from a observational data set. A possible 
alternation is to adopt the weight factor: 

w(s,k) = ; ^— — — - - . (31) 

This modified weight factor does not minimize the errors, however, it might be useful in a 
realistic situation of data analysis. The weight factor contains (Vo(k)), then it must be solved 
with equation (15) as a coupled system. However, we can easily solve the coupled equation 
numerically by iterations. In this case the amplitude of the covariance matrix (AVo(k) 2 ) 
is different from the expression (30). It is obtained by substituting (31) into (24), but the 
expression is rather complicated than (30). As we show below, the use of the modified 
weight factor alters the expression of the Fisher matrix element. However, the difference is 
not significant. For example, the use of the modified weight factor increases the 1 — o error 
of the equation of state Aw (see below) by 10 % level. Then our conclusion does not depend 
on the choice of the above weight factors. 



2.2. Fisher Matrix Element 

In order to estimate the accuracy to which we can constrain cosmological parameters 
with a measurement of the power spectrum, we employ the Fisher matrix approach. With 
the Fisher-matrix analysis, one can estimate the best statistical errors on parameters from 
a given data set (Kendall & Stuart 1969). For this reason, this approach is widely used to 
estimate how accurately cosmological parameters are determined from large surveys, such as 
the large scale structure of galaxies, cosmic microwave background anisotropy measurement 
and supernova data sets. (See e.g. Tegmark, Taylor, & Heavens 1997, Jungman et al. 1996a, 
1996b, Zaroubi et al. 1995, Fisher et al. 1995, Tegmark et al. 1998b). In general, the Fisher 
matrix is defined by 

where L is the probability distribution function of a data set, given model parameters 6*j. For 
a measurement of galaxy power spectrum, the Fisher matrix is presented using a suitable 
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approximation in the references (Tegmark 1997, Tegmark et al. 1998a). Here we briefly 
review a derivation of the Fisher matrix: For simplicity, we adopt the approximation of the 
Gaussian probability distribution function for A"P(k), 



L oc exp 



\J dk J dk'AV(k)C(k,k')AV(k') 



(33) 



where C(k, k') is the inverse matrix of the covariance matrix (A'P(k)A'P(k')). Then, using 
(25), we have 



where n(k) is 



u K min 1 J 



_ fdsn(s)Ms,kr(P (k,s) + l/n(s)) 2 
K{ } ~ [J rfsn(s)2^(s,A;)2]2 ' [6b) 



and (Vo(k)) is the expectation value of the power spectrum, i.e. expression (15). The Fisher 
matrix depends on the weight factor, and we have 

-^ M / fe lu^Mf (36) 

for the optimal weight factor (29), and 

An[JdssMsni + n(s)(Vo(k)))-r (r . 
{) /rf SS 2 n(s) 4 (P (A;,s) + l/n(s)) 2 (l + n( s )(Po(A;)))- 4 1 ' 

for the modified weight factor (31), respectively, where we assumed that the mean number 
density is a function of the distance s(— |s|), and Ail is the survey area. The result is con- 
sistent with the previous work: In the limit that the linear redshift distortion, the geometric 
distortion, and the light-cone effect are switched off, the above results reproduce the previous 
work (e.g., Tegmark 1997). By using the Bayes theorem, the probability distribution in the 
parameter space can be written 

pie,) oc exp Y,{0i - 0?)fm - ef)] , (38) 

ij 

where we have assumed that errors in the target model parameters 8f are small. Thus 
the Fisher matrix gives the uncertainties in the parameter spaces, which are described by a 
Gaussian distribution function around Of '. 
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3. 



CONSTRAINTS 



In this section we investigate a prospect of the power spectrum analysis of the LRG 
sample in the SDSS. As pointed out by Matsubara & Szalay (2002), the LRG sample can 
be a useful tool to constrain the cosmological parameters. Here we assess the potential of 
the LRG sample to constrain the equation of state of the dark energy. Because the LRG 
sample is distributed out to a redshift ~ 0.5, the geometric distortion is substantial. This is 
the reason why the power spectrum analysis of the LRG sample can constrain the equation 
of state wx even if the original matter power spectrum (or the transfer function) does not 
depend on wx- In comparison to the LRG sample, the SDSS quasars are distributed out to 
z ~ 5, however, the spatial distribution is very sparse. Then the constraint from the quasar 
sample will not be very tight due to the large shot-noise contribution (Matsubara & Szalay 
2002, Yamamoto 2002). 

The clustering properties of the LRG sample have not been analyzed so far, then we 
here assume simple linear bias models for the clustering. (See equation (46) and section 3.3 
for a luminosity dependent bias model). To evaluate the Fisher matrix element, the number 
density of the sample n(z) is an important factor. In the present paper, we adopt the LRG 
sample with the comoving number density in the reference by Eisenstein et al. (2001, Figure 
12 in their paper), as a function of the redshift in the range 0.2 < z < 0.55. The typical 
number density of LRGs is n ~ 10 -4 [/i^Mpc] -3 . The peak value of the power spectrum is 
10 5 [/i^Mpc] 3 . Thus n(z)P (k,z) does not exceed 10. 



We restrict ourselves to a spatially flat universe for modeling the cosmology. This 
assumption can be justified by the inflationary universe scenario and the recent results of 
the cosmic microwave anisotropy measurements. We consider a cosmological model with the 
dark energy component with the equation of state which is variable in time. We assume that 
the time variation of the equation of state wx(z) = Px/px is slow. The equation of motion, 
d(p x V(z)) + p x dV(z) = with a volume V(z) oc (1 + z)~ 3 , yields 



3.1. Modeling the Power Spectrum 




(39) 



With the parameters w and u, we assume wx can be parametrized by 




(40) 
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where z* is a constant, then we have 

/*(*) = (1 + -) 3 -p[( ITZF (^ )J- (41) 

Denoting the matter density parameter by fl m , the comoving distance is 

r{z) = -L [ Z = (42) 

H o Jo y/n m (l + z') 3 + (l-n m )f x (z') 

We can model the power spectrum in redshift space s(z) as (Ballinger, Peacock, & 
Heavens 1996) 



P„(M*)) = I d^( q , - ^ - ^^'4 ( 43 ) 

where P ga i(q\\, q±, z) is the galaxy power spectrum, q n and q ± are the wave number components 
parallel and perpendicular to the line-of-sight direction in real space with r(z), and we define 

C ' (Z) = M7) and c ^ ) = 4Ty (44) 

The geometric distortion is described by scaling the wave number by c^(z) and c ± (z). Then, 
even if the original matter power spectrum or the transfer function does not depend on 
the parameter of the dark energy wx, the power spectrum observed in redshift-space can 
depend on w X - Thus we can test the nature of dark energy. This idea traces back to the 
geometric test pointed out by Alcock & Paczynski (1979). We model the power spectrum in 
the distribution of galaxy by 

Pgnfa,^,*) = (l + ^^)V) 2 P mass (g,z), (45) 

where q 2 = q 2 + q1, b(z) is a scale independent bias factor, P mass (q, z) is the CDM mass power 
spectrum, and we defined f(z) = d\nD(z)/d\na(z) with the linear growth rate D(z) and 
the scale factor a(z). The term in proportion to f(z) describes the linear distortion (Kaiser 
1987). Following the work by Matsubara & Szalay (2002), we work within the linear theory 
of density perturbation, adopting the fitting formula of the transfer function by Eisenstein 
& Hu (1998). Here we consider a simple bias model parametrized by, 

6 (,) = l + ( foo -l)_^, (46) 

where bo and p are constants. For f(z) we adopt the fitting formula developed by Wang and 
Steinhardt (1998) 

f(z) ~ n(z) a{wx) , (47) 



- 12 - 



with 



a(w x ) 



3 3 (l-w x )(l-3w x /2) 



(i-n(z)), 



(48) 



5-w x /(l-w x ) 125 (l-6w x /5) 3 



fi m + (l-fi m )(l + z)-3/ x (z)' 



(49) 



In the present paper, we work within the linear theory of density perturbations and we do 
not consider the nonlinear and the finger of Got effects. The nonlinear effect is substantial on 
small scales in the nonlinear regime (Mo, Jing & Boerner 1997; Magira, Jing & Suto 2000). 
The inclusion of the nonlinear effect will not alter our result significantly, however, any 
uncertainty of modeling the nonlinear effect reduces the capability to determine cosmological 
parameters precisely (cf. Watkins et al. 2002). 



In our analysis we focus on the sensitivity of the power spectrum analysis on the nature 
of the dark energy, then we consider the parameters (fi m , w, v, bo,p)- We assume that other 
cosmological parameters such as the baryon density fib, c"8 and the index of the initial density 
power spectrum n are well determined. We might expect to obtain such information in the 
near future from other observations of the cosmic microwave background anisotropies and 
the large scale structure of the main galaxy sample in SDSS, etc. Actually the WMAP 
result have demonstrated the cosmological parameters can be determined from the cosmic 
microwave background anisotropies. In the present paper, we take Vt m = 0.28, Q b h 2 = 0.024, 
h = 0.7, cr s = 0.9 and n = 1 (Spergel et al. 2003). On the other hand, the bias is an 
annoying problem, though the physical mechanism has been extensively investigated (e.g. 
Mo & White 1996). Recent investigation reports that the bias depends on the luminosity 
and galaxy type (Norberg et al. 2001, 2002). Thus the bias depends on sample, and it should 
be determined from the same data simultaneously. In our investigation, we here consider 
the marginalized probability function integrating over the bias parameters bo and p in (38). 

Figure 1 shows the contour of the marginalized probability function on fl m — w space, 
which is obtained by integrating (38) with respect to (u, bo, p). Here, we adopted the 
target parameters (Of) as follows: Q m = 0.28, w = —1, v = 0, b = 1.5 and p — 1. The 
target parameter of the equation of state corresponds to the cosmological constant. The 
dependence of the bias modeling is discussed in more detail in section 3.3. In (40), we here 
set = 0.13. We also assume the survey area 10 4 deg 2 , which corresponds to the complete 



3.2. Constraints on wx 
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SDSS sample. 2 We adopted the range of the integration in (34) being k min = 0.001 /iMpc -1 
and /c max = 1 ftMpc -1 . Our results do not depend on k min as long as k min £ 0.01 /iMpc -1 , 
while depending slightly on fc max , but our conclusions are not significantly altered as long as 
fc max ^ 0.3 /iMpc -1 . It is clear from Figure 1 that the power spectrum is more sensitive to 
Q m than w. In our modeling the transfer function depends on Q m , but not on wx- The effect 
of wx on the power spectrum comes from the geometric distortion. The weaker dependence 
of the power spectrum on w than Vt m originates from this fact. 

Figure 2 shows the contour of the marginalized probability function on w — v space 
obtained by integrating (38) with respect to (Q m , bo, p). The target parameters are same 
as those in Figure 1. Similarly, we set = 0.13 in (40). In this case, it is clear from Figure 
2 that the degeneracy between w and v is broken. Note that w = wx(z*). Thus the power 
spectrum analysis of the LRG sample gives the equation of state around z = 0.13, almost 
independently from v. Conversely, the LRG sample is not very sensitive to probe v. 

The minimum error in determining the equation of state is Aw ~ 0.1, which is given 
by integrating the probability function over the other parameters (fl m , v, b , p). The error 
can be reduced, when Q m is determined from other observation independently. The result 
depends on the target parameters. Figure 3 shows the I — a error Aw as function of fi^j. The 
solid curves assume w tT = —1, while the dashed curves assume w tr = —0.8. The other target 
parameters are same as those in Figure 1. The error is sensitive to the density parameter 
Q^, however, is not sensitive to w tr . The sensitivity to the equation of state deteriorates in 
higher matter density universes. This is because the fraction of the dark energy to the total 
energy decreases as the matter density becomes higher. As long as £ 0.3, however, the 
minimum error is Aw ^,0.1. 



Here we consider a more realistic bias model. Recent investigations with the 2df galaxy 
redshift survey report that the clustering amplitude depends on galaxy type and luminosity 
(Norberg et al. 2001, 2002). The selection effect from an apparent limiting magnitude can 
have significant influences on the redshift evolution of the bias. Denoting the luminosity 
dependent bias by b(z,L), we define the luminosity-averaged bias as 



3.3. A Luminosity Dependent Bias Model 



L 



(X! 



(p( y L)b( y z,L)dL 



b(z) 




<p(L)dL 



(50) 



2 The dependence on the survey area will be discussed in the final section. 
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where (f>(L) is the luminosity function. In modeling b(z,L), we assume that b(z,L) can be 
approximately written as b(z,L) = b(z)b(z = 0,L). Namely, the luminosity dependence of 
the bias does not depend on the redshift. Here we also assume b(z) to be described by (cf. 
Mo & White 1996) 

which is normalized to yield 1 at z — 0, i.e. b(0) = 1. For b(z = 0, L), we adopt the result 
by Norberg et al. (2002), who found that the luminosity dependence can be fitted by 

b(z = 0,L) X A + (1-A) — , (52) 

where A is constant. Norberg et al. report A = 0.8 for early-type galaxy (2002). Combining 
(51) and (52), we assume b(z,L) can be written in the form 

^■ L ) = ( 1 + ^)( A+ < 1 -^)- (53) 

Then, expression (50) yields 

where F((3,x) is the incomplete Gamma function and x(z) = Li im (z)/L*, and we assumed 
the Schechter luminosity function 

«^ = ^£)*«p(- (55) 

In the present paper, we adopt the fitting formula for the luminosity function in the reference 
by Madgwick et al. (2002) for early-type galaxies (a = —0.54 and — 5 log 10 h = —19.58), 
and m = 19.5 for the apparent limiting magnitude. Figure 4 shows b(z) as a function of z, 
where we set bo = 1.2 and A = 0.8, and we have not considered the luminosity evolution and 
the k-correction for simplicity. 

We consider the bias parameterized by (54) with b and A, instead of (46) with b 
and p, and repeat the evaluation of the Fisher matrix. Here we fix the target parameters, 
Q m = 0.28, w — — 1 and v — 0. Figure 5 shows the marginalized best statistical error Aw, 
integrating over the parameters except for w as function of the target parameters bo with 
A fixed A = 0.6 (short-dashed curve), A = 0.8 (solid curve), A = 0.9 (long-dashed curve) 
and A = 1 (dot-dashed curve). Thus, when the bias evolution can be fitted in the form of 
(54), Aw ~ 0.1, which is almost the same level as that in Figure 3 for A <, 0.9. However, 
for the case A — 1, the error Aw increases out to 0.15, depending on b . This suggests that 
modeling of the redshift-evolution of the bias can be problematic to constrain the equation 
of state of dark energy using the LRGs. 
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4. OPTIMIZED SAMPLE 

In this section, we address the problem of determining the optimized sample in order 
to best constrain the equation of state. This problem might be only of theoretical interest, 
but having such information could be useful in planning a survey of galaxies or clusters of 
galaxies. The same problem has been considered with regard to supernova data and lensing 
systems by several authors (Huterer & Turner 2001, Spergel & Starkman 2002, Yamamoto 
et al. 2001). To address the problem stated above, we introduce the Fisher matrix element 
per object as follows: When a survey volume is small enough so that the light-cone (redshift 
evolution) effect is negligible, (34) is written 

1 /" fcmax AVn 2 dP (k,s)dP (k,s) 3 



F " - L (i+«p.(M))» «r — aer k dla *• (56) 

where we assumed n to be constant, and AV is the survey volume. Then Fij/AV can be 
regarded as the Fisher matrix element per unit volume, and similarly Fij/(AVn) is regarded 
as the Fisher matrix element per object. In a strict sense we must choose k m \ n and /c max 
depending on AV and n. However, we here fix k m - m = 0.01 /iMpc -1 and /c max — 1 ftMpc -1 . 
Then we investigate Fij/(AVn) as a function of z and n. Here note that s is the function of 
the redshift, i.e. s = s(z). 

In this section, for simplicity, we consider the model in which the equation of state is 
constant wx(z) = w, which is equivalent to the assumption fx(z) — (1 + z) 3( - 1+w \ Regarding 
the Fisher matrix as the 2x2 matrix corresponding to the elements w and b , we consider 
the marginalized Fisher matrix element integrating the probability function over bg'. 

F- 

Fww F ww — . (57) 

Figure 6 shows [AVn,/ F^] 1 / 2 as a function of z with n fixed asn= 10~ 2 /i 3 Mpc -3 (dashed 
curve), 10~ 4 ^, 3 Mpc _3 (solid curve), and 10~ 6 /i 3 Mpc -3 (long dashed curve). Figure 6 in- 
dicates that objects in the range 0.4 ^ z £ 1 are most efficient for measuring w, because 
[AVn/ F^yj] 1 / 2 is regarded as the error in determining w per object at the redshift z. Note 



that the error using the total objects Aw scales as 1/yN = 1/yfiAV. Figure 7 shows 
[AVn/ F^w] 1 / 2 as a function of n with the redshift fixed at z — 0.6. This figure indicates that 
the number density around n = 10^ 4 /i 3 Mpc~ 3 is most efficient for measuring w. For very 
sparse samples, like quasars, which have a typical number density n = 10~ 6 /i 3 Mpc~ 3 , the 
errors in measuring parameters using the power spectrum is large due to the shot-noise. On 
the other hand, for a sample with high number density, the shot-noise is of minor importance, 
but the efficiency in constraining w per object decreases as well. Thus the objects in the 
range of redshift 0.4 £ z £ 1 with the number density n ~ 10~ 4 /i 3 Mpc~ 3 are the optimized 
sample. 
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5. CONCLUSIONS 

In summary we derived a rigorous optimal weighting scheme for power spectrum analy- 
sis, which is useful for samples in which the light-cone effect and the redshift-space distortions 
are substantial. Our result is a simple extension of the work by FKP, and we obtained a 
generalized optimal weight factor which minimizes the errors of the power spectrum estima- 
tor. As an application of our formula, we investigated the capability of the LRG sample in 
SDSS to constrain the equation of state parameters by evaluating Fisher matrix elements. 
Even if the transfer function for the matter power spectrum does not depend on wx, the 
power spectrum analysis of the redshift-space sample can constrain wx due to the geomet- 
ric distortion. To incorporate uncertainties of redshift evolution of the clustering bias, we 
considered the marginalized probability function by integrating over the parameters of the 
bias models. This analysis shows that the LRG sample in SDSS has a serviceable potential 
for constraining the equation of state around z = 0.13 with 1 sigma errors at the 10 % 
level, if other fundamental parameters are well determined in an independent fashion. We 
also showed that this conclusion is not altered in the case of a bias model incorporating the 
redshift-evolution due to selection effect depending on luminosity. However, even in our real- 
istic treatment of the bias, we made simplifications: Uncertainties including the stochasticity 
and the nonlinearity in modeling the bias are not considered in our investigation. Then, tests 
on the bias properties will be required for more definite conclusions. 

In the present paper, we assumed 10 4 deg 2 as the complete SDSS survey area. When 
the planed survey area are not achieved, the capability to constrain the equation of state 
reduces. As the Fisher matrix element is in proportion to the survey area AQ, then the 
statistical error Aw increases in proportion to VAf2. For example, when we assume 5 x 10 3 
deg 2 and 7.5 x 10 3 deg 2 as the final SDSS survey area, Aw increases by 30 % and 15 %, 
respectively, as long as the inhomogeneity of the incomplete survey area does not cause 
additional systematic errors. 

In section 4, we considered the optimized sample to constrain wx using power spectrum 
analysis. We found that it is most advantageous to have the sample with the comoving 
number density n ~ 10~ 4 /i 3 Mpc~ 3 in the range of redshift 0.4 £ z £ I. For such a sample, 
the efficiency per object to constrain wx is optimized. Information from anisotropic power 
spectrum would improve the capability of constraining the parameters, as demonstrated in 
the 2df QSO sample (Outram et al. 2001). 
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Fig. 1. — Contour of the marginalized probability function in Q m -w space. The curves show 
1 — a, 2 — a and 3 — a contours. The target parameters are Q m = 0.28, w = —1, v = 0, 
b = 1.5 and p — 1. In (40) we set = 0.13. The survey area 10 4 deg 2 is assumed. 
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Fig. 2. — Contour of the marginalized probability function in w-v space. The meaning of 
the curves and the parameters are same as those in Figure 1. 
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Fig. 3. — The statistical 1 — a error Aw as a function of the target parameter Q^. The solid 
and dashed curves adopt w tr = — 1 and w tr = —0.8, respectively. For each pair of curves, the 
lower curve assumes that Q m is determined by other independent method, while the upper 
curve assumes that fl m is determined simultaneously by the power spectrum analysis. The 
other target parameters are same as those in Figure 1. 
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Fig. 4. — b(z) as a function of the redshift. We adopted b — 1.2 and A = 0.8. 
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Fig. 5. — The statistical error Aw as a function of the target parameter 6 with fixing 
A = 0.6 (short-dashed curve), A = 0.8 (solid curve), A = 0.9 (long-dashed curve) and A = 1 
(dot-dashed curve) in the bias model in section 3.3. The other target parameters are same 
as those in Figure 1. Here we assume that fl m is determined simultaneously by the power 
spectrum analysis. 
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Fig. 6. — [AVn/ Fww] 1 ^ 2 as function of z with n fixed. The dashed curve, the solid curve, 
and the long dashed curve, assume n = 1CT 2 /i 3 Mpc~ 3 , 1CT 4 /i 3 Mpc -3 , and 1CT 6 /i 3 Mpc~ 3 , 
respectively. The target parameters are same as those in Figure 1. 
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Fig. 7. — [AVn/Fww] 1 / 2 as a function of n. Here the redshift is fixed at z — 0.6. The target 
parameters are same as those in Figure 1. 



